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32 ■ Abstract 



We investigate the power-suppressed corrections to fragmentation functions in 
flavour singlet deep inelastic lepton scattering, to complement the previous 
results for the non-singlet contribution. Our method is a dispersive approach 
based on an analysis of Feynman graphs containing massive gluons. As in 
non-singlet deep inelastic scattering we find that the leading corrections are 
proportional to l/Q 2 - 
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1 Introduction 



The study of fragmentation functions in deep inelastic lepton scattering (DIS) has received 
a great impetus from the increasing quantity and kinematic range of the HERA data 
IJ. While these functions cannot be calculated perturbatively in QCD, their asymptotic 
scaling violations (logarithmic Q 2 dependence) can be calculated, enabling measurements 
of the strong coupling constant a s . In addition there are power-suppressed (higher-twist) 
contributions to the scaling dependence, which need to be estimated in order to make use 
of the wide Q 2 coverage of HERA. 

Recently so-called 'renormalon' or 'dispersive' methods of estimating power-suppressed 
terms have been formulated: for descriptions see p, [3|, |fl. Such approaches have been 
applied with some success to e + e~ fragmentation functions || and event shape variables 
||, and to flavour non-singlet DIS structure functions 0, [7], §, fragmentation functions || 
and event shape variables JHJ. Note that special care has to be taken when considering 



non-inclusive observables such as fragmentation functions and event shape variables [11 



There are also renormalon model results for photon-photon scattering [|12| and structure 



functions in flavour singlet DIS [13 



In the present paper we extend the dispersive method to fragmentation functions in 



flavour singlet DIS, using many of the same techniques as in [13]. As in non-singlet DIS 
and e + e~ annihilation, the predicted leading power corrections to these quantities are 
proportional to 1/Q 2 , but their functional forms are different. In particular there are 
contributions proportional to (logQ 2 )/Q 2 , which are not found in the non-singlet case. 
While the magnitude of the correction is not known, the hypothesis that it is related to a 
universal low-energy strong coupling implies that it is given by universal non-perturbative 
parameters. 



In the following section we review the approach of references |3], |13j . Section |3] presents 



the standard leading-order perturbative treatment of DIS fragmentation functions. In 
section [| we estimate the power-suppressed corrections using the method outlined in section 
^|. Our results are summarized briefly in section [5]. 



2 The Dispersive Approach to Power Corrections 

We assume that the QCD running coupling a s (k 2 ) can be defined for all positive k 2 , and 
that apart from a branch cut along the negative real axis there are no singularities in the 
complex plane. It follows that we may write the formal dispersion relation: 

a s (k 2 ) = - jT '-^Pst/i 2 ), (2.1) 
where the 'spectral function' p s represents the discontinuity across the cut: 

Ps(/i 2 ) = -^(a s (/iV") - a s (p 2 e- m )\ = -^Disc « s (-/i 2 ) • (2.2) 

Non-perturbative effects at long distances are expected to give rise to a non-perturbative 
modification to the perturbatively-calculated strong coupling at low scales, <5a s (/i 2 ) = 
a s (/x 2 ) — af T (/x 2 ), a^ T (/i 2 ) being the perturbatively-calculated running coupling [|J. This 
modification gives rise to power-behaved corrections to QCD observables. (Note that here 
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T (/z 2 ) refers only to the next-to-leading order perturbatively-calculated running coupling, 
and so is itself well-behaved down to low scales, without any divergences: the Landau pole 
appears when we include an arbitrary number of loop insertions in the propagator. Hence 
5a s ([i 2 ) is a well-defined quantity for all positive /i 2 .) 

We now consider the calculation of some observable F in an improved one-loop approx- 
imation which takes into account one-gluon contributions plus those higher-order terms 
that lead to the running of a s . We obtain the characteristic function T{p? : /Q 2 ), which is a 
one-loop evaluation of F but with the gluon mass set equal to \i. We may then write |§ 



f°^W) \H» 2 /Q 2 )-m] . (2.3) 

Jo fi Z 1 



The behaviour of T{e) at small e generates a correction to F, which we denote SF. Terms 
of J-'ie) which are analytic at e = do not contribute to 5F but non-analytic terms at 
small e do contribute, the relevant terms being: 

F-a^yfe =► ^ = --7T> (2.4) 

and ^ 

^ ~ a2 2n el ° ge ^ 5F = a2 ~<^' 

where 



A q = ^l ^4^5a s (fi 2 ) . (2.6) 



— / — p, q da s (^ ' 
lix Jo p z 

Studies of the non-singlet contribution to DIS structure functions suggest that Ai 
0.2GeV 2 §. 

This is the formulation of the dispersive approach to power behaved corrections used 
where there is a single gluon propagator. However, the calculation which follows in section 
f| involves two gluons, and so the above argument needs to be generalised. 

Each gluon has an associated dispersive variable, so we obtain a characteristic function 
in both of these, jF(ei,e 2 ). Following the above argument, we require those terms of T 
which are non-analytic at in both arguments. However, in the case where there are two 
internal gluons constrained to have the same 4-momentum k, we can simplify this to require 
only one dispersive variable. By defining p = —k 2 /Q 2 , we see that the dependence of the 
characteristic function T on e\ and e 2 is given by 

H^)=j , + dP (!- P \ (2.7) 
J (P + ei)(p + e 2 ) 

where the integration limits and the function / depend on the particular calculation. This 
may be expressed in partial fractions in the form 

= fl iM^M, (2 . 8) 

ei - e 2 

where 
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Now the generalisation of the one-gluon result ( |2.3|) to the multiple-gluon case (which 
can easily be seen from a generalisation of the argument in ||) is 

^^4 Ps ^l)p s ^l)[^\/Q\ - T^\IQ\ 0) - .F(0, fij/Q 2 ) + .F(0, 0)] . 

Q H X fi 2 

(2.10) 

Substituting fl2.8|) into this gives, as shown in |13J, that the relevant non-analytic terms in 
jF(e) are: 

D 1 

T ~ aieloge =^ 5F = a x — , 



and 

T ~ \a 2 t log 2 e 
where D\ and D 2 are defined by: 



^ = a 2 -log- 



Dl = 



dp 2 

o /i 2 



p 2 (2a,Gu 2 )oa s (/i 2 ) - [oa s (/i 2 )] ; 



lo gj D 2 = [ < ^-p 2 \ogp 2 (2a s (p 2 )Sa s (p 2 ) - [5a s (p 2 )f 

U\ JO [L v 



(2.11) 
(2.12) 

(2.13) 
(2.14) 



While we expect the form of a s (fj, 2 ), and hence D\ and D 2 , to be universal, we have as 
yet no numerical values for them, (unlike for A2, which is defined by (2J3) with q = 2). It 
will be necessary therefore to extract values for D\ and D 2 , either from experimental results 
or from some model of the form of a s (p 2 ) (of which various models have been proposed 

I)- 



3 DIS Fragmentation Functions 

We consider the deep inelastic scattering of a lepton of momentum / from a nucleon of 
momentum P, with transfer of a virtual photon of momentum q. The standard variables 
for describing DIS are Q = \J—q 2 (Q > 0), the Bjorken variable x = Q 2 /2P-q (0 < x < 1), 
and y = P ■ q/P ■ I ~ Q 2 /xs, where s is the total cm. energy squared, (0 < y < 1). 
We wish to study the fragmentation function F h (z), defined by 

{Z ' ,X,Q ' = dxdQ 2 dz I dx^Q 2 ' ' ( } 

for a given hadron species h, as a function of the variable z = 2ph ■ q/q 2 , where pt is 
the momentum of the resultant hadron. Since the fragmentation products of the remnant 
of the nucleon are expected to be travelling in directions close to that of the incoming 
nucleon, i.e. in the 'remnant hemisphere' p z < 0, we consider only hadrons produced in 
the 'current hemisphere' p z > 0. Such hadrons are expected to be fragmentation products 
of the scattered parton. Thus z takes values < z < 1. 



The differential cross sections in ( |3.1|) can be decomposed into 
d 3 a h 2na 2 



{[1 + (1 - y) 2 }F£(x, z) + 2(1 - y)F£(x, z)} (3.2) 
dxdQ 2 Q4 {t 1 + ( X ~ v) 2 ]Mx) + 2(1 - y)F L (x)} , (3.3) 



dxdQ 2 dz Q 4 
d 2 a 2ira 2 
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Figure 1: Flavour singlet contribution to deep inelastic scattering 



where Ft and Fl are the transverse and longitudinal structure functions and Fj> and 
F£ are their generalisations to include fragmentation. They are related to F\ and F<i by 
F T {x) = 2F 1 (x) and F L (x) = F 2 (x)/x - 2F 1 (x). 

The parton model assumptions are that p = xP/C, (x < £ < 1) and ph = zp p /C, 
{z < ( < 1), where p p is the momentum of the fragmenting parton (i.e. for the diagram 
in figure [I], p p = r for quark fragmentation and p p = k for antiquark fragmentation). Thus 
we get 

F?(x,z) = £ X: /f £f3WZ)DZ(z/QF iJk (Z,C) (3-4) 



incoming outgoing 
partons j partons fc 



X 



E tjj^imM) (3.5) 



incoming 
partons j 



where i = T or L, j(x/£) is the parton distribution function for the parton j, D^(z/() is 
the fragmentation function for the parton k into the hadron h, and i^jjfc^, C) an d 
are the (generalised) parton level structure functions. Note that all the functions in the 
above equations have a scaling violation or Q 2 dependence which is not shown explicitly. 
(For simplicity we are neglecting any contribution from weak interactions, i.e. Z° or W 1 * 1 
exchange.) 

In the parton model, to order a°, we have 

F*(x,z) = E e2 M^)D h q (z)+q(x)D h M (3.6) 
q 

F h L {x,z) = (3.7) 

Ft{x) = J2e 2 q [q(x) + q(x)} (3.8) 
q 

F L {x) = (3.9) 
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where q(x) and q(x) are the quark and antiquark distributions in the target nucleon. 

There are 0(a s ) contributions arising both from real and virtual gluon emission, as 
described in ||. Here we concentrate on the singlet contribution, i.e. the contribution 
from the gluon content of the target hadron, as shown in figure [I]. There are no virtual 
corrections to the singlet contribution at this order. 



4 Power Corrections in Flavour Singlet DIS 

In the normal perturbative treatment of DIS, the asymptotic freedom of QCD enables 
us to treat the initial state partons as free particles confined within the nucleon; and so 
in a singlet calculation we would start from a free gluon and convolute the perturbative 
result with the gluon distribution function g(x). We do not know how to do this in a 
calculation of power corrections, since the models we use consider modifications to the 
gluon propagator (loop insertions in the renormalon model, or, equivalently, a 'mass' in the 
dispersive approach). Let us therefore perform the calculation by considering our initial 
state gluon to be radiated from a fermionic parton. We may then try to recover the singlet 
contribution to the power corrections by deconvoluting the result with the quark to gluon 
splitting function, as performed in |14| , or we may leave the result as it is and interpret it 
as a second-order non-singlet contribution. We might hope that these two interpretations 
would give similar predictions for power-suppressed corrections. 




Figure 2: Diagram generating flavour singlet contribution 



Consider the contribution to DIS from the diagrams in figure 0. It is convenient to work 



in the Breit frame of reference ||15|| , which is the rest-frame of 2xP + q. In this frame the 
momentum transfer q is purely spacelike, and we choose to align it along the +z axis. 

The momentum of the initial state parton isp = xP/£, (x < £ < 1); and let us introduce 
the variables p = —s 2 /Q 2 , rj = P ■ r / P ■ q, fj = P ■ kj P ■ q , x the azimuthal angle between 
r and s, and 7 the azimuthal angle between r and k. There is also an overall azimuthal 
angle <p. 
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In the Breit frame the kinematics are given by: 

(4.1) 
(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 

The definitions of p, rj and fj along with the on-shell conditions for the outgoing particles 
require that 



p = 


§Q(l/z,0,0,-l/z) 


p = 


|g(i/e,o, o,-i/o 


Q = 


§Q(0,0,0,2) 


s = 


±<2(s ,s_l cosx, s±sinx, s 3 ) 


r = 


±Q(zo,z±,0,z 3 ) 


k = 


\Q{z G , z± cos 7, z± sin 7, z 3 ) 



1 > 4 

zo=r]+z i 



1 , s i 



Conservation of the 0th and 3rd components of 4-momentum give the conditions 



V + V + pi 

z\ z\ s\ 

— + — H — 

4i] 4r] 4p£ 



(4.7) 
(4.8) 
(4.9) 

(4.10) 
(4.11) 



while conservation of transverse momentum requires that s±, z± and z± satisfy the triangle 
inequalities 



\z± — s±\ < z± 
\z± — z±\ < s± 



(4.12) 
(4.13) 



Two variables, a and f3, are required to parametrise the permitted values of s±, z±_ and 
z±. Let us choose to write: 



4^(1-^(1-0 
z 1 = 



4^(1-0 
4p(l - 

1-77 



[1 — a)pt; + ar/rj — 2 cos/3v/a(l — a)r]rjp^ 



(1 — a)rj + ar/p^ + 2cos/3v/o;(l — a)r]r]p^ 



(4.14) 
(4.15) 
(4.16) 



We could equally have used a and /3, defined similarly but with the outgoing quark and 
antiquark interchanged. In particular, if z\ = 4afj(l — 77) (1 — 0/£> then 



a = 



(l-rj)(l-fj) 



(1 — a)p£ + arji] — 2cosf3da(l — a)r]r]p^ 



(4.17) 



Given s±, z± and z±, the angles x an d 7 are determined up to a sign. We may then 
choose a, ft, 77 and p as the independent variables, with phase space 



< a < 1, < (3 < 71, < 77 < 1, < p < (1 - 77)/^. 



(4.18) 
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To estimate power corrections to perturbative calculations, we must perform the calcu- 
lations as though the gluons had small non-zero masses p\ = EiQ 2 and /i 2 , = e 2 Q 2 . Using 
the machinery of section |2| we may write 

w x e 1 J r i(x,z;e 1 )-e 2 J r i(x,z;e 2 ) . 
Fi(x, z\ ei, e 2 ) = , (4.19) 

£l — ^2 

where, by the symmetry between the outgoing quark and antiquark, 



H*> *\ <0 = E E 4 [ f jf f c; eM*/o [i#(*/c) + ^/o; , 

(4.20) 



9 q' 



and C is given in terms of the other variables by 



2r-q _ £ + q(l-0 a(l - Q 

? 2 e 17 e 



c = fL^ = iXif^-I^ - ^ViZ^ . (4.21) 



The quantities Cj are the integrated perturbative matrix elements with modified gluon 
propagators. 

To integrate the matrix elements we apply the operator 

Integrating away p' and making substitutions for r and k gives 

Q 2 f z ± dz ± d(pdr] z ± dz ± d~fdfj 5 ( A + V B + 2z ± z ± cos7 y 

, (4.zo) 



32(27r) 5 7 r) fj y/B + 2z ± z ± cos-f 

where A and B do not depend on 7. 

Next we integrate over 7. There are two values satisfying the delta-function condition, 
and they differ by a sign. This gives 



Q 2 f drj dr] dz^dz 2 



32(27r) 5 J r] 7] 2z±z±\ sin7| 
Applying the parametrisation in terms of a and f3, we find that 

d(zlzl) 



(4.24) 



n r] = 32(1 -0 2 sin PyJa(l-aWf}3p/? (4.25) 
22 ± 2 ± |sin7| = 8(1 - £) an/3y/a(l - at)ijfjp/Z . (4.26) 
Therefore the integral operator is 

Q 2 (i-0 



8(2tt) 



J dpdr]dadf3d(J) ■ (4.27) 



We are interested to find the differential coefficient functions as a function of (. In the 
current hemisphere (( > 0) there is a singularity at £ = 1, associated with the antiquark 
becoming soft or collinear with the proton. This divergence must be investigated separately. 
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The divergences associated with the quark becoming soft or collinear occur when ( < 0, 
and so are not seen in the current hemisphere. 
For < £ < 1, then we have 



Git, C; e) = --^S^ 1 f da d P r d(3 r n ' ( 428 ) 

e 8(27r) & jo Jo Jo Jo £ + a(l — £,) p + e 

where the function /(£, a, (3, (f), p) is obtained from the appropriate matrix elements, 
excluding the denominators of the gluon propagators. We calculate and Cl using 
projections of the hadronic tensor representing the QCD part of the process. Since we wrote 
the kinematics without ^-dependence, the projections obtain ^-dependence to compensate. 
We find that the integration over is trivial for C L and elementary for C T - 

The part of this that is non-analytic as e — > is the contribution to the integral near 
p = 0. Therefore we proceed by expressing / as a series expansion in pa about 0, up 
to 0(p 2 ). The integral over (3 is then straightforward, being a sum of terms of the form 
const, cos™ f3, and upon performing this integral the terms in non-integer powers of p vanish. 
Then note that the non-analytic parts of the integrals over p are given by: 



n 

1 dp -> (-ly-Vloge, (4.29) 



o p + e 
p n log p 



o p + e 



dp -> \{-l) n - 1 e n \og 2 e. (4.30) 



The integral over a is straightforward. 

To study the divergence at ( = 1, consider, for some small A and some arbitrary 
function <?(C), 



A-A 



dCCite C; e)0(C) 



e 8(2tt)5 Jo Jo h J^-l^ 'Jo Jo > + 



In order to produce valid series expansions in p, it transpires that we have to divide the 
i] integral into two regions: 

W 1 ~ g+i(i-0 - V ~ 1 ~ p ^ ~ p / K ' and 
(ii) 1 - p£ - p/« < T) < 1 - p^, 

where /t is some small arbitrary quantity on which the final answer should not depend. (In 
practice we perform the calculation in the limit of small k.) 

For region (i) we may directly expand / in p without any problems; and since A is small 
we can expand g(() as 

9(0 = 9(1) ~ $ + Q ^~ (l - VW(1) + ^(A 2 ) . (4.32) 

Note that because of the limits we have chosen for the r\ integral, contributions up to 
0(p 2 ) will also arise from all higher terms in the series expansion of /. These contributions 
vanish as k — > 0, so the calculation performed in this limit is valid. 
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For region (ii) such a naive power series becomes invalid. We use the symmetry between 
the outgoing quark and antiquark to see that /(£, 77, a, (3, 0, p) = /(£, fj, a, /3, 0, p) and that 
the integration measure drjdad/3 = dfjdad/3. We may therefore write the required integral 
as: 



lQ 2 (l-0 ti 



da 



dW d/?/ d<j>^—g((). 
Jo Jo p + e 



e 8(2vr) 5 

We may now write 7] = Xp and then expand / as before. We may also expand g as: 

£ + fi(i-0 



0(0=0(1: 



;i-^y(i)+o(A 2 



(4.33) 



(4.34) 



In contrast to the previous situations there are now no terms in non-integer powers of p. 
The integration over /3 may be performed by writing t = tan(/3/2), and that over a by 
making a series of substitutions. The integration over A (i.e. if) can then be performed, 
and, since we are only concerned about the limit k — > 0, we neglect all terms that vanish 
in this limit. The p integration then proceeds as above. 

Having calculated the contributions near £ = 1 that do not vanish as A — > 0, we may 
express the divergence in terms of delta functions and plus prescriptions by comparison 
with the following results: 



0(0 



l-A(l-C)- 

9(0 

1-A (1 " 01 



-d( 
-d( 



<7(1) logA + 0(A) 
9(1) 



(ftl) log A + (/(!) + 0(A) 



1-A 
1 

1-A 



g{Q6{l-QdC 
g(()6'(l - ()d( 



9(1) 
9(1) 



(4.35) 
(4.36) 
(4.37) 
(4.38) 



Putting all this together we have 



Ct(£,0 



rf(£, C) + G^(0S(i - + h$ } (0 ( |*(i - log(e0 



(0), 



(i-O- 



loge 



+ r«(e, + 4 1] (0^(i - + 4 1} (0 (Va - iog(eo - 7^ 



+4 2) (0^'(l - + #H0 ( 1^(1 - log(eO + 



r(2) 



4(i-p 2 (i + 2e 2 

r(l) 



loge + 



+#no 1^(1-0 io g (eo 



(i-Oi 

r^^O + c? (0^(1-0 
1 



o+; 

eloge (4.39) 



(1-0- 



e log e , 



(4.40) 



where the functions F do not diverge at £ = 1. The explicit forms of the functions T, G 
and i7 are given in the appendix. 

The terms that diverge as e — > are responsible for the logarithmic scaling violations 
to the structure function (c.f. ||), while the terms proportional to eloge and elog 2 e give 
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the 1/Q 2 power corrections. We then obtain power corrections to the generalised structure 
functions, given by 



where 



5F l (x) 



c 



Q 2 (2vr) 2 J x i 



(4.41) 



4 2) (0 



+ 



(i-C) + (l-C) 



+ 



g£> (0+^(0 log 
4 2) (0 + 4 2) (0iog 



Q 2 
Q 2 



(DUz) + zDiJ(z)) (4.42) 



c 



^(0 

:i-o+ 



+ 



4 1} (0 + ^ ; (0iog 



£fcw(*/C) 



(4.43) 



and the charge-weighted sum over all quark and antiquark fragmentation functions is 

DiM = J2e 2 q [D h q (z) + D h q (z)}, (4.44) 



while the total quark distribution in the target hadron is 



(4.45) 



For simplicity, let us suppose that D q (z) is the same for all quarks and antiquarks. This 
approximation is valid if there is negligible heavy flavour production and if we sum over an 
appropriate set of final-state hadrons. Using the definition of the fragmentation function 



JJ) and using the order a° results (|3.6|) - ( |3.9|) , we find that 
q c . w (x)5F h (z,x;Q 2 ) = (6F*(x, z)-D h q (z)6F T (x))+ ^'^^ F^x, z)-D h q (z)5F L (x)) 

where the charge-weighted quark distribution in the target hadron is 

Qc-w(x) = X)eJ[g(aO + q(x)} . 



(4.46) 
(4.47) 



as 



The power corrections to the totally inclusive structure functions are given in [|T^, [14 

D h q (z)5 Fi ( X ) = ^l^^/ 1 f toWQSCti) , (4.48; 
where 

D 2 £ 



5C T (£) = 2C? (0 + 2^(0 log 
6C L (0 = 2Ci 1} (0 + 2^(0 k)g 



eQ 2 
eQ 2 



(4.49) 
(4.50) 
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This leads to a power correction to the fragmentation function given by 



r ZTi/i/ ^\ D 1 T R C F 

6F {z ^ Q) = q^W? 



(4.51) 



where the transverse and longitudinal contributions are 



Kn 



1 dj q tot (x/Cj 



+ ( G?(0 - C^(0 + ^(0(1 - 1 logO ) DUz 



c 

.(1) 



4 2) (0 + 4 2) (0 log ^ } (D h c _ w (z) + zD h J( 



Q 2 



(4.52) 



Kj 



* £ <7c-w i x ) 



i;fh^-^)n L{ z /0 



+ ( 4 1} (0 - (£) + 4^ (0(i - 1 logO ) 



0- 



(4.53) 



and the explicit forms of the coefficient functions T, C, G and H may be found in the 
appendix. The quantities D\ and D 2 are defined in equations ( |2.13| ) and ( |2.14j ). Note that 
while there are both 1/Q 2 and (logQ 2 )/Q 2 elements in the transverse contribution, the 
(logQ 2 )/Q 2 terms cancel in the longitudinal contribution. 



5 Results and Conclusions 

The 1/Q 2 power corrections arising from the diagrams shown in figure ||] are given above 
in equations ( |4.51| ), ( [4.52[ ) and ( [4.53| ) . While we do not know the values of D 1 and D 2 , 



we can still make some qualitative predictions. Figure |3] shows plots of Kt(z) and K^{z) 
at fixed x-values of 0.001, 0.01, 0.1 and 0.5. A logarithmic scale is used for the z axis. 
These were calculated at Q 2 = 500 GeV 2 , using the corresponding MRST (central gluon) 
parton distributions ||16|| , and ALEPH parametrisations for the light quark fragmentation 
functions [O. The value of D 2 was set to be 0.06 GeV 2 , i.e. approximately A 2 , following the 



approach of [13| and [Q. (The qualitative behaviour of the Ki does not change provided 
we keep D 2 <tiQ 2 .) 

The parton distributions and fragmentation functions used in the above calculation also 
have weak Q 2 dependences, behaving as logQ 2 , which are not shown explicitly in equa- 
tions ( [4.52j ) and ( |4.53| ) but which combine with the explicit power corrections to produce 
additional scaling violations behaving as (logQ 2 )/Q 2 - They are implicit within the distri- 
bution or fragmentation function and any calculation of power corrections should take the 
scale-evolution of these functions into account. 

As can be seen from the graphs, Kt{z) and Kl(z) both tend to zero at large z, and 
both are divergent and negative at small z, with Kt dominating over Kl- The magnitudes 
of Kt and are larger at smaller x, although the qualitative behaviour does not change 
much with x provided x is kept reasonably small. Therefore if Di is positive, we have a 
negative correction to F h at small z. If Di is negative then the correction is positive. 
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Figure 3: Graphs showing Kt(z) and Kl(z) for various fixed values of x. 



We do not know the sign or magnitude of Di, and there is no a priori reason why it 
should be either positive or negative. However, since A2 is positive, we presumably have 
from (2.6) that 5a s (fi 2 ) is predominantly positive. Considering the definition Q2.13 ) of D\, 



and noting that 2a s 5a s — (8a s ) 2 = 2a^ T 8a s + (5a s ) 2 , we therefore might conjecture that 
D\ would also be positive. Preliminary results from HERA may suggest a large negative 



correction at small z, lending support to this hypothesis [18 
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Appendix 

Those functions which for brevity were not written in full in the main body of this paper 
are given below. Note that the functions T do not diverge as £ — > 1. 



4 0) (0C) 



+ 



4(1 + 2(1 



0i-0 4 e(i-C) 3 £(i-C) 



2 , 2(1 + Q , 



C 



log 



+ 



4(1-0 2(1 



e 



ea-0 3 ea-o 2 301-0 



+ 



4-8^ + 8Q-4Q-3Q log^ 
3£ 



(A.l) 



12 



r?(OC) 



180 



oi-o 6 e(i-o 5 

4 + 19£ + 4p 2£(l + £)" 



192(1 + + 6 ( 9 + 34 £ + 9 ^ 2 ) 4 ( x + 14 £ + 14 ^ 2 + £ 3 



(i-C) 2 



1-C 



log 



oi-o 4 



oi-o 3 

180(1-0 102(1-0) 



18(1 + 3£ - 30 - 0) 2(2 - 7£ + 70 - 20 



Oi -C) 3 Oi -C) 2 

20(1-2^ + 20) 2(2-2^ + 0) 



Oi -C) 5 Oi -C) 4 

(1-0(4 + 70 + 40) 

«i-C) 



if (00 



i-e+ec 

144 



oi -o 6 ea-c) 



ec 

192(1 + 8(7 + 30^ + 90) 



(A.2) 



Oi -C) 4 



8(8 + 10£ + Q 

(i-C) 3 



4(5£ + 20 



40 



(i-O 2 i-C 

8(1 + 11^-90-3^ 



log 



+ 



' i-e+ec ^j 

g(l + e-2Q 

(i-C) 2 



144(1-0 120(1-0) 



+ 



Oi -C) 5 Oi -C) 4 

4(2-150 + 100 + 30)" 

1501 -0 



(i 



c 



(1 



G 



(o 



G 



^(2 



G 



H. 



(o 



(i 



(2 



(1 



2 + 250 - 250 - 20 + 150 log f + 150 log f 

5£ 

4(17 + 750 - 1250 + 330 + 30 log£ + 750 log £ - 450 logf) 

225£ 

26 - 90£ + 90Q - 26Q + 12 log g - 9Qo g e ~ 9plog£ + 12g 3 logg 

6 + 20£ + 900 - 900 - 200 - 60 + 4 50 log £ + 750 log f 

7 - 20£ + 20Q - 7Q - 24Q log£ 

12£ 

-2(92 + 1500 - 4 25p + 1830 + 30 log£ - 225p log£ 
+ 1500 log f + 450 log / (2250 

4 + 3£ - 3Q - 4Q + 6£ log £ + 6Q log £ 
3£ 

-05-5£ + 21og£ + 2£logp) 

eioge 

4(2 - 150 + 10O + 30 - 150 logO 
15i ' 



(A.3) 

(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 

(A.9) 

(A.10) 
(A.ll) 
(A.12) 

(A.13) 



References 



[1] ZEUS Collaboration, M. Derrick et al.: Z. Phys. C 67 (1995) 93, ( |hep-cx 9501Trq ). 
HI Collaboration, S. Aid et al: Nucl. Phys. B 445 (1995) 3, ( |hep-ex 9505003|) . 
C. Adloffet al: Nucl. Phys. B 504 (1997) 3, ( |hep-ex 9707005] ) . 

[2] For reviews and classic references see V.I. Zakharov: Nucl. Phys. B 385 (1992) 452, 
and A.H. Mueller, in QCD 20 Years Later, vol. 1 (World Scientific, Singapore, 1993). 



13 



M. Beneke, V.M. Braun and V.I. Zakharov: Phys. Rev. Lett. 73 (1994) 3058, ( |hep-ph 
94053041 ) 



P. Ball, M. Beneke and V.M. Braun: Nucl. Phys. B 452 (1995) 563, ( [hep-ph 95023001) . 
M. Beneke and V.M. Braun: Phys. Lett. B 348 (1995) 513, ( frep-ph 94112291) . 
M. Beneke and V.M. Braun: Nu cl. Phys. B 454 (1995) 253, ([hep-ph 9506452j ). 
M. Beneke: CERN-TH/98-233, (|hep-ph 98074431) . 



M. Neubert: Phys. Rev. D 51 (1995) 5924, (|hep-ph 9412"265| ). 



Yu.L. Dokshitzer and N.G. Uraltsev: Phys. Lett. B 380 (1996) 141, ( [hep-ph 9512401) . 
G. Grunberg: CPTH-PC463-0896, (|hep-ph 9608375] ). 



G. Grunberg: Ecole Polytechnique preprint CPTH-S505-0597, ( hep-ph 9705"290[) . 

[3] Yu.L. Dokshitzer, G. Marchesini and B.R. Webber: Nucl. Phys. B 469 (1996) 93, 
Qhep-ph 95123361) . 



[4] Yu.L. Dokshitzer, V.A. Khoze and S.I. Troyan: Phys. Rev. D 53 (1996) 89, ( |hep-ph 
95064251 ). 



D.V. Shirkov and I.L. Solovstov: Phys. Rev. Lett. 79 (1997) 1209, ( [hep-ph 9704331) 
G. Grunberg: Phys. Lett. B 372 (1996) 121, (|hep-ph 9512"203l) 



D.V. Shirkov and I.L. Solovtsov: JINR Rapid Comm. 2 [76] (1996) 5, ([hep-ph 



96043631 ). 

D.V. Shirkov: Nucl. Phys. Proc. Suppl. 64 (1998) 106, ( hep-ph 9708480|) . 
A.I. Alekseev and B.A. Arbuzov: Mod. Phys. Lett. A 13 (1998) 1747, ( |hep-ph 



97042281 ). 



B.R. Webber: JEEP 10 (1998) 012, ( |hep-ph 9805181 ). 



[5] M. Dasgupta and B.R. Webber: Nucl. Phys. B 484 (1997) 247, (|hep-ph 9608394|) . 



P. Nason and B.R. Webber: Phys. Lett. B 395 (1997) 355, flhep-ph 9612353; ). 
M. Beneke, V.M. Braun and L. Magnea, Nucl. Phys. B 497 (1997) 297, ( |hep-ph 
97013091 ). 



[6] B.R. Webber: Phys. Lett. B 339 (1994) 147. Yu.L. Dokshitzer and B.R. Webber: 
Phys. Lett. B 352 (1995) 451, ( [hep-ph 95042191) . 

G.P Korchemsky and G. Sterman: Nucl. Phys. B 437 (1995) 415, ( |hep-ph 941l2TT| ). 



G.P. Korchemsky and G. Sterman: in the Proceedings of the 30th Rencontres de 

Moriond, Meribel les Allues, France, March 1995, ( |hep-ph 950539 1| ). 

G.P. Korchemsky, G. Oderda and G. Sterman: to appear in the Proceedings of the 

5th International Workshop on DIS and QCD, ([hep-ph 9708346|) . 

R. Akhoury and V.I. Zakharov: Phys. Lett. B 357 (1995) 646, ( |hep-ph 95042i8| ). 



R. Akhoury and V.I. Zakharov: Nucl. Phys. B 465 (1996) 295, ( hep-ph 9507253]) 



P. Nason and M.H. Seymour: Nucl. Phys. B 454 (1995) 291, ( |hep-ph 95063T7D 



Yu.L. Dokshitzer and B.R. Webber: Phys. Lett. B 404 (1997) 321, (|hep-ph 9704298Q . 



[7] M. Dasgupta and B.R. Webber: Phys. Lett. B 382 (1996) 273, ( [hep-ph 96043gg| ). 

[8] E. Stein, M. Meyer-Hermann, L. Mankiewicz and A. Schafer: Phys. Lett. B 376 (1996) 
177, ( |hep-ph 960135B1) . 



M. Meyer-Hermann, M. Maul, L. Mankiewicz, E. Stein and A. Schafer: Phys. Lett. B 
383 (1996) 463, ( [hep-ph 9605229)) ; ibid. 393 (1997) 487 (E). 



M. Maul, E. Stein, A. Schafer and L. Mankiewicz, Phys. Lett. B 401 (1997) 100, 



14 



( Ihep-ph 96123001) . 

M. Maul, E. Stein, L. Mankiewicz, M. Meyer- Hermann and A. Schafer, ( |hep-ph 
9710392| ). 



M. Meyer-Hermann and A. Schafer, ( [hep-ph 9709349| ). 

[9] M. Dasgupta, G.E. Smye and B.R. Webber: JEEP 04 (1998) 017, (|hep-ph 9803382| ). 

[10] M. Dasgupta and B.R. Webber: Eur. Phys. J. CI (1998) 539, (|hep-ph 97042971) . 

M. Dasgupta and B.R. Webber: Cambridge preprint Cavendish-HEP-98/02, ( |hep-ph 



98092471 ). 



[11] Yu.L. Dokshitzer, A. Lucenti, G. Marchesini and G.P. Salam: Nucl. Phys. B 511 
(1998) 396, ( [hep-ph 9707352]) . 

Yu.L. Dokshitzer, A. Lucenti, G. Marchesini and G.P. Salam: JEEP 05 (1998) 003, 
( Ihep-ph 98023811) . 



G.P. Salam: Milan preprint IFUM-623-FT, ( |hep-ph 9805^|) . 
[12] F. Hautmann: Phys. Rev. Lett. 80 (1998) 3198, ( |hep-ph 97102561) . 
[13] G.E. Smye: Cavendish-HEP-98/14, (|hep-ph 9810292ft . 

[14] E. Stein, M. Maul, L. Mankiewicz and A. Schafer: Torino preprint DFTT 13/97, to 
appear in Nucl. Phys. B , (|hep-ph 9803342|) . 

[15] R.D. Peccei and R. Riickl: Nucl. Phys. B 162 (1980) 125. 

K.H. Streng, T.F. Walsh and P.M. Zerwas: Z. Phys. C 2 (1979) 237. 

L.V. Gribov, Yu.L. Dokshitzer, S.I. Troyan and V.A. Khoze: Sov. Phys. JETP 68 

(1988) 1303. 

[16] A.D. Martin, R.G. Roberts, W.J. Stirling and R.S. Thorne: Eur. Phys. J. CI (1998) 
463, ( |hep-ph 9803l45|) . 

[17] ALEPH Collaboration, D. Buskulic et al: Phys. Lett. B 357 (1995) 487, Phys. Lett. B 
364 (1995) 247 (E). 

[18] HI Collaboration: submission 531 to the 29th International Conference on High- 
Energy Physics, Vancouver, 1998. 

ZEUS Collaboration: submission 809 to the 29th International Conference on High- 
Energy Physics, Vancouver, 1998. 



15 



